Fuel limiting method in diesel engines having exhaust gas recirculation

ABSTRACT

A method of controlling the fueling rate of a compression ignition engine having an exhaust gas recirculation system (EGR) and a turbocharger. The method includes the steps of predicting the exhaust gas air/fuel ratio as a function engine operating characteristics and the operator-requested fueling rate. This value is compared to an air/fuel limit value stored in a table of values indexed by engine operating characteristics. The air/fuel limit is the value below which visible smoke occurs. The predicted air/fuel ratio is compared to the air/fuel ratio limit and either (1) the fuel is delivered at the requested rate, (2) the start of injection timing is modified and the fuel is delivered at the requested rate or (3) the fueling rate is limited as a function of the requested fueling rate and measured or determined engine operating parameters.

TECHNICAL FIELD

This invention relates to compression ignition engines having exhaust gas recirculation systems and, more particularly, to methods of limiting fuel delivery to avoid visible smoke in turbocharged diesel engines equipped with exhaust gas recirculation (EGR) systems.

BACKGROUND OF THE INVENTION

High performance, high speed diesel engines are often equipped with turbochargers to increase power density over a wider engine operating range, and EGR systems to reduce the production of NOx emissions.

Turbochargers use a portion of the exhaust gas energy to increase the mass of the air charge delivered to the engine combustion chambers. The larger mass of air can be burned with a larger quantity of fuel, thereby resulting in increased power and torque as compared to naturally aspirated engines.

A typical turbocharger consists of a compressor and turbine coupled by a common shaft. The exhaust gas drives the turbine which drives the compressor which, in turn, compresses ambient air and directs it into the intake manifold. Variable geometry turbochargers (VGT) allow the intake airflow to be optimized over a range of engine speeds. This is accomplished by changing the angle of the inlet guide vanes on the turbine stator. An optimal position for the inlet guide vanes is determined from a combination of desired torque response, fuel economy, and emissions requirements.

EGR systems are used to reduce NOx emissions by increasing the dilution fraction in the intake manifold. EGR is typically accomplished with an EGR valve that connects the intake manifold and the exhaust manifold. In the cylinders, the recirculated exhaust gas acts as an inert gas, thus lowering the flame and in-cylinder gas temperature and, hence, decreasing the formation of NOx. On the other hand, the recirculated exhaust gas displaces fresh air and reduces the air-to-fuel ratio of the in-cylinder mixture.

Both the VGT and EGR regulate gas flow from the exhaust manifold, and their effect is, therefore, coupled through the conditions in the exhaust manifold. Excessive EGR rates displace the intake of fresh air and may lead to incomplete combustion of the injected fuel which, in turn, could cause visible levels of smoke and increased levels of emissions. Moreover, this could negatively affect fuel economy and/or performance. Thus, for effective control of diesel engines with EGR systems, it is necessary to control the EGR flow precisely, not only in steady state but also in transient conditions. Precise EGR control requires monitoring the EGR flow, the fresh air flow, and the intake mixture composition to control the combustion process and, thereby, avoid smoke production and particulate emissions.

Conventionally, fuel demand is calculated as a function of engine speed and accelerator pedal position, and fueling rate limiters restrict the actual amount of fuel injected as a calibrated function of intake manifold pressure (MAP), intake manifold temperature, and engine speed (RPM) to avoid smoke. Thus, such systems operate without regard to measurements or estimates of the amount of fresh air or oxygen available in the intake manifold. In particular, fuel limiting is conventionally accomplished as a function of charge density, without regard to charge composition, or as a function of compressor mass airflow, without regard to the additional oxygen contributed by the EGR flow.

Knowledge of the amount of fresh air available for combustion in the intake manifold, however, is important. In a diesel engine, the generated torque is directly proportional to the amount of injected fuel, provided enough air is available. Visible smoke is also strongly related to the air/fuel ratio. Consequently, it is desirable to deliver fuel to the engine at a rate which generates the torque demanded by the driver, yet maintains the air/fuel ratio above the threshold at which visible smoke occurs.

DISCLOSURE OF THE INVENTION

One object of the present invention is to provide a fast airflow response while maintaining the air/fuel ratio above the threshold at which smoke occurs.

Another object is to control fuel injection quantity and timing to avoid smoke generation.

The invention is advantageous in that it improves system performance by dynamically estimating intake manifold oxygen concentration, allowing for faster airflow response without visible smoke and hence, reduced turbo lag.

According to the present invention, the foregoing and other objects and advantages are obtained by a method of controlling the fueling rate of a compression ignition engine having an EGR system and a turbocharger. The method comprises the steps of generating a turbocharger intake airflow value (W_(c1) (k)), an EGR flow value (W₂₁ (k)) and engine intake airflow value (W_(1e) (k)). From these measured or estimated values, the intake burnt gas fraction is calculated as a function of W_(c1) (k) and W₂₁ (k). Once the intake burnt gas fraction value is obtained, the start of fuel injection signal and requested fueling rate are determined as a function of pedal position, engine speed and possibly other engine operating conditions. A predicted exhaust gas air/fuel ratio (λ_(e) (k)) is then generated as a function of W_(1e) (k) and the requested fueling rate. Furthermore, an air/fuel ratio limit λ_(lim) (k) is determined from a stored lookup table based on engine operating conditions. This is the air/fuel ratio limit below which visible exhaust smoke will occur. The predicted exhaust gas air/fuel ratio is compared to the exhaust gas air/fuel ratio limit (λ_(lim) (k)) necessary to avoid smoke. If λ_(e) (k)≧λ_(lim) (k) fuel is delivered at the requested rate. Otherwise, the start of fuel injection value is modified to increase the exhaust gas air/fuel ratio limit. This new exhaust gas air/fuel ratio limit is then compared to λ_(e) (k). If λ_(e) (k)≧λ_(lim),new (k), the start of injection timing value is modified and fuel is delivered at the requested rate, otherwise, the fuel is limited as a function of the intake burnt gas fraction, W_(ie) (k) and λ_(lim),new (k).

Other objects and advantages of the invention will become apparent upon reading the following detailed description and appended claims, and upon reference to the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the invention, reference should now be had to the embodiments illustrated in greater detail in the accompanying drawings and described below by way of examples of the invention. In the drawings:

FIG. 1 is a schematic view of a compression ignition engine system having an EGR system and a VGT in accordance with one embodiment of the present invention;

FIG. 2 is a logic diagram describing a method of controlling the engine fueling rate in accordance with one embodiment of the present invention;

FIGS. 3a-c are graphs illustrating estimated values according to one embodiment of the fuel limiting method; and

FIGS. 4a-e are graphs illustrating the performance of the fuel limiting method of one embodiment of the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENT

Turning first to FIG. 1, there is shown a simplified schematic diagram of a compression ignition engine system 10 equipped with an exhaust gas recirculation (EGR) system 12 and a variable geometry turbocharger (VGT) 14. A representative engine block 16 is shown having four combustion chambers 18. Each of the combustion chambers 18 includes a direct-injection fuel injector 20. The duty cycle of the fuel injectors 20 is determined by the engine control unit (ECU) 24 and transmitted along signal line 22. Air enters the combustion chambers 18 through the intake manifold 26, and combustion gases are exhausted through the exhaust manifold 28 in the direction of arrow 30.

To reduce the level of NOx emissions, the engine is equipped with an EGR system 12. The EGR system 12 comprises a conduit 32 connecting the exhaust manifold 28 to the intake manifold 26. This allows a portion of the exhaust gases to be circulated from the exhaust manifold 28 to the inake manifold 26 in the direction of arrow 31. An EGR valve 34 regulates the amount of exhaust gas recirculated from the exhaust manifold 28. In the combustion chambers, the recirculated exhaust gas acts as an inert gas, thus lowering the flame and in-cylinder gas temperature and decreasing the formation of NOx. On the other hand, the recirculated exhaust gas displaces fresh air and reduces the air-to-fuel ratio of the in-cylinder mixture.

The turbocharger 14 uses exhaust gas energy to increase the mass of the aircharge delivered to the engine combustion chambers 18. The exhaust gas flowing in the direction of arrow 30 drives the turbocharger 14. This larger mass of air can be burned with a larger quantity of fuel, resulting in more torque and power as compared to naturally aspirated, non-turbocharged engines.

The turbocharger 14 consists of a compressor 36 and a turbine 38 coupled by a common shaft 40. The exhaust gas 30 drives the turbine 38 which drives the compressor 36 which, in turn, compresses ambient air 42 and directs it (arrow 43) into the intake manifold 26. The VGT 14 can be modified during engine operation by varying the turbine flow area and the angle at which the exhaust gas 30 is directed at the turbine blades. This is accomplished by changing the angle of the inlet guide vanes 44 on the turbine 38. The optimal position for the inlet guide vanes 44 is determined from the desired engine operating characteristics at various engine speeds.

As can be appreciated from FIG. 1, both the EGR 12 and the VGT 14 regulate gas flow from the exhaust manifold 28. The effect of the EGR and VGT is, therefore, jointly dependent upon the conditions in the exhaust manifold 28.

All of the engine systems, including the EGR 12, VGT 14 and fuel injectors 20 are controlled by the ECU. For example, signal 46 from the ECU 24 regulates the EGR valve position, and signal 48 regulates the position of the VGT guide vanes 44.

In the ECU 24, the command signals 46, 48 to the EGR 12 and VGT 14 actuators are calculated from measured variables and engine operating parameters by means of a control algorithm. Sensors and calibratable lookup tables provide the ECU 24 with engine operating information. For example, manifold absolute pressure (MAP) sensor 50 provides a signal 52 to the ECU 24 indicative of the pressure in the intake manifold 26. Likewise, exhaust manifold pressure (EXMP) sensor 54 provides an EXMP signal 56 to the ECU 24 indicative of the pressure in the exhaust manifold 28. Further, an aircharge temperature sensor 58 provides a signal 60 to the ECU 24 indicative of the temperature of the intake aircharge 42. A mass airflow (MAF) sensor 64 also provides a signal 66 indicative of the compressor intake airflow to the ECU 24. Additional sensory inputs can also be received by the ECU along signal line 62 such as engine coolant temperature, engine speed, and throttle position. Additional operator inputs 68 are received along signal 70 such as acceleration pedal position. Based on the sensory inputs and engine mapping data stored in memory, the ECU controls the EGR valve to regulate the intake airflow (MAF), controls the VGT to regulate the intake manifold pressure (MAP) and controls injectors 20 to regulate fuel delivery.

FIG. 2 describes the logic flow of the ECU to regulate fuel delivery to avoid the production of smoke.

Throughout the specification, the following notations are used in describing measured or calculated variables:

    ______________________________________                                          p                pressure (kPa)                                               T                       temperature (K)                                        m                       mass (kg)                                              W                       mass flow (kg/s)                                       F                       burnt gas fraction                                     α.sub.egr     EGR valve position                                         ______________________________________                                    

Furthermore, the following subscripts are used to denote regions of the engine system:

    ______________________________________                                         1                intake manifold 26                                            2                   exhaust manifold 28                                        e                   engine block 16                                            t                   turbine 38                                                 c                   compressor 36                                              ______________________________________                                    

Finally, the following thermodynamic constants are referenced for air at 300 K: ##EQU1## Hence, the symbol p₁, for example, refers to the intake manifold pressure in kPa. Combined subscripts, such as "e2", refer to flows from the first to the second subsystem.

FIG. 2 describes the logic routine to accomplish fuel limiting based on estimates of intake and exhaust aircharge composition. This logic routine resides in the ECU memory and is executed as part of the routine used to control the timing and duty cycle of the fuel injectors 20.

In step 102, the VGT compressor airflow (W_(c1) (k)), EGR flow rate (W₂₁ (k)) and engine intake flow (W_(1e) (k)) are computed.

The compressor airflow (W_(c1) (k)) can be measured by MAF sensor 64 (FIG. 1) , or can be estimated based on measurements of intake manifold pressure from MAP sensor 50, exhaust manifold pressure from EXMP sensor 54, and the EGR valve position.

The EGR flow value (W₂₁ (k)) is calculated as a function of intake manifold pressure (p₁), exhaust manifold pressure (p₂), the EGR temperature (T₂₁), and EGR valve position (α_(egr))according to the following equation:

    W.sub.21 (k)=ƒ.sub.1 (α.sub.egr)p.sub.2 /(RT.sub.21).sup.1/2 φ(p.sub.1 /p.sub.2)                                   (1)

The EGR temperature (T₂₁) is determined from a steady-state map based on engine operating conditions. Alternatively, T₂₁ can be assumed to be a constant. The air charge temperature (T_(c1)) can also be measured by a temperature sensor such as sensor 58 of FIG. 1, or estimated based on engine operating conditions.

In equation (1), f₁ (α_(egr)) represents the effective flow area of the EGR valve as a function of the position of the EGR valve, R represents the difference between the pressure specific heat constant and volume specific heat constant, and φ represents a standard orifice equation having the following form:

    φ(r)=[(2γ/(γ-1))(r.sup.2/γ -r.sup.(γ+1)/γ ].sup.1/2  for r≧(2/(γ+1)).sup.γ/(γ-1)

    φ(r)=γ.sup.1/2 (2/(γ+1)).sup.(γ+1)/(2(γ-1)) for r<(2/(γ+1)).sup.γ/(γ-1)

The flow out of the intake manifold into the combustion chambers (W_(1e) (k)) is then given by:

    W.sub.1e (k)=(NV.sub.d n.sub.vol /(RT.sub.1 120))p.sub.1   (2)

wherein N represents the engine speed, V_(d) represents the displacement volume of the engine, and n_(vol) represents the volumetric efficiency of the engine. The volumetric efficiency is stored in the ECU memory as a function of one or more of the following variables: intake manifold pressure, intake manifold temperature, fuel rate, engine speed, and engine coolant temperature.

Given a measured value of the compressor flow rate (W_(c1) (k)), the EGR flow rate (W₂₁ (k)) can also be estimated. The time rate of change of the intake pressure can be expressed as:

    dp.sub.1 /dt=-Ap.sub.1 +(R.sub.γ /V.sub.1)(W.sub.c1 T.sub.c1 +W.sub.21 T.sub.21)                                       (3)

where

    A=(NV.sub.d γn.sub.vol /(V.sub.1 120))               (4)

Thus, from equation (3), the only unknown term is the EGR flow rate W₂₁. An observer is, therefore, constructed to dynamically estimate the product W₂₁ T₂₁ by interpreting the scaled enthalpy flow as the state of a dynamic system whose dynamics are assumed to be zero. Assume the state of the estimator is the scaled enthalpy flow, z=W₂₁ T₂₁, and letting p₁ and z be estimates of intake manifold pressure and z, respectively, a pressure error term and flow error term can be defined as follows:

    e.sub.p =p.sub.1 -p.sub.1

    e.sub.z =z-z

The following differential equations are then observers for manifold pressure and scaled enthalpy flow:

    dp.sub.1 /dt=-(A+MA)p.sub.1 +MAp.sub.1 +(R.sub.65 /V.sub.1)(W.sub.c1 T.sub.c1 +z)                                              (5)

    dz/dt=L(p.sub.1 -p.sub.1)                                  (6)

where L and M are calibratable constants whose value is >0. The error dynamics for equations (5) and (6) are:

    de.sub.p /dt=-(1+M)Ae.sub.p +(R.sub.65 /V.sub.1)(e.sub.z)  (7)

    de.sub.z /dt=-Le.sub.p                                     (8)

Appropriate values for the design parameters M and L can be determined; a requirement being that the eigenvalues of the error system are in the left half complex plane. For example, assuming M=0.5, L=0.5, and V₁ =0.003 m³, and the engine operating conditions are yielding n_(vol) =0.85, the eigenvalues are given by -7.3353 and -4.5647.

The estimate of the EGR flow value (W₂₁) is then given by the following equation:

    W.sub.21 =z/T.sub.21                                       (9)

In order to implement equation (9) in the digital ECU, it can be discretized with a sufficiently small sampling period δt. In such a case, the value of W₂₁ is governed by the following equations:

    p.sub.1 (k+1)=p.sub.1 (k)+δt[-(A+MA)p.sub.1 (k)+MAp.sub.1 (k)+(R.sub.γ /V.sub.1) . . . ##EQU2##

Alternatively, the EGR flow value can be defined on the basis of the ideal gas law--instead of the first law of thermodynamics as in equations (4) and (12)--as follows: ##EQU3## where M' and L' may be different from M and L.

In accordance with another embodiment of the invention, W₂₁ can be calculated by another method based on the first law of thermodynamics wherein the following equation defines the intake manifold pressure:

    dp.sub.1 /dt=(R.sub.γ /V.sub.1)(W.sub.c1 T.sub.c1 +W.sub.21 T.sub.21 -W.sub.1e T.sub.1)                                        (16)

Applying a Laplace transform to both sides of equation (16) and multiplying equation (16) by 1/(s/τ+1) results in the following equation:

    sp.sub.1 /(s/τ+1)=(R.sub.γ /V.sub.1)(1/(s/τ+1))(W.sub.c1 T.sub.c1 +W.sub.21 T.sub.21 -W.sub.1e T.sub.1)            (17)

From equation (17) estimates for the time rate of change of the intake pressure and EGR flow can be defined as follows:

    dp.sub.1 /dt=sp.sub.1 /(s/τ+1)                         (18)

    W.sub.21 =W.sub.21 /(s/τ+1)                            (19)

    f=(1/(s/τ+1))(W.sub.c1 T.sub.c1 -W.sub.1e T.sub.1)     (20)

Substituting these values in equation (17), the filtered EGR flow, W₂₁, is defined as:

    W.sub.21 =(1/T.sub.21)((V.sub.1 /(γR))dp.sub.1 /dt-f))(21)

Intake pressure, p₁, and intake aircharge temperature, T_(c1), are measured values from MAP sensor 50 and temperature sensor 58 of FIG. 1. The remaining variables are known or can be resolved. For example, the intake mass airflow, W_(c1), is obtained from MAF sensor 64. Similarly, the engine intake flow rate, W_(1e), is obtained from the mapped volumetric efficiency, measured intake manifold pressure, and engine speed as in equation (2). Also, the EGR temperature T₂₁, can be taken as a constant, or mapped as a function of measured engine operating conditions. Finally, the intake manifold temperature, T₁, is obtained from the steady state equation:

    T.sub.1 =(W.sub.c1 T.sub.c1 +W.sub.21 T.sub.21)/(W.sub.c1 +W.sub.21)(22)

Preferably, to implement the control logic in the digital ECU, the logic can be sampled over discrete time periods, δt, resulting in the following controller equations:

    d(k+1)=d(k)+δt(-τd(k)+p.sub.1 (k))               (23)

    dp.sub.1 (k)/dt=τ(p.sub.1 (k)-τd(k))               (24)

    f(k+1)=f(k)+δt(-τf(k)+τ(W.sub.c1 (k)T.sub.c1 (k)-W.sub.1e (k)T.sub.1 (k))                                           (25)

    W.sub.21 (k)=(1/T.sub.21 (k))((V.sub.1 /(γR))(dp.sub.1 (k)/dt)-f(k))(26)

wherein V₁ represents the volume of the intake manifold, dp₁ (k)/dt represents an estimate of the time rate of change of the intake manifold pressure, f(k) represents the filtered difference between the compressor mass air flow enthalphy and engine intake flow enthalpy, and T₂₁ (k) represents the EGR temperature.

As a further embodiment, the EGR flow value can be defined on the basis of the ideal gas law, instead of the first law of thermodynamics as in equations (16), (25) and (26), as follows:

    d(k+1)=d(k)+δt(-τd(k)+p.sub.1 (k))               (27)

    dp.sub.1 (k)/dt=τ(p.sub.1 (k)-τd(k))               (28)

    f(k+1)=f(k)+δt(-τf(k)+τ(W.sub.c1 (k)-W.sub.1e (k)))(29)

    W.sub.21 (k)=(V.sub.1 /RT.sub.1 (k)) dp.sub.1 (k)/dt-f(k)  (30)

wherein V₁ represents the volume of the intake manifold, T₁ (k) represents the temperature of the intake manifold, dp₁ (k)/dt represents an estimate of the time rate of change of the intake manifold pressure, and f(k) represents the filtered difference between the compressor mass flow rate and engine intake flow rate.

The performance of the EGR flow rate estimator as defined by equation (12) is illustrated in FIG. 3b. FIG. 3b shows a graph of estimated W₂₁ using equation (12) (line 202) versus the simulated actual W₂₁ (line 200) over a period of 200 seconds.

From step 102, it is advantageous to have knowledge of the EGR mass flow rate W₂₁ because it allows for a more precise, or less conservative, fuel limiting scheme. Specifically, the composition of the intake aircharge and the amount of oxygen contributed by the EGR flow are taken into account in deciding whether to limit the fueling rate in order to avoid the production of smoke. This process is carried out by the ECU in steps 104 through 122.

At step 104, the intake burnt gas fraction F₁ is estimated. The differential equation governing the dynamics of the burnt gas fraction in the intake manifold is based on a mass balance and is defined as:

    dF.sub.1 /dt=1/m.sub.1 [-(W.sub.c1 +W.sub.21)F.sub.1 +W.sub.21 F.sub.2 ](31)

The mass in the intake manifold is obtained using the ideal gas law, and substituting in equation (31) to obtain:

    dF.sub.1 /dt=(RT.sub.1 /(p.sub.1 V.sub.1))[-(W.sub.c1 +W.sub.21)F.sub.1 +W.sub.21 F.sub.2 ]                                       (32)

It is evident that equation (32) is stable. Hence, a direct integration of the right hand side yields a stable open-loop observer for F₁. The results of this observer for the intake burnt gas fraction are demonstrated in FIG. 3c. In FIG. 3c, line 208 graphs the estimated intake burnt gas fraction from equation (32), and line 210 represents the simulated actual intake burnt gas fraction.

In step 106, the start of fuel injection (SOI(k)) value is determined as a function of engine speed and fueling rate from a lookup table stored in the ECU memory.

In step 108, the requested fueling rate (W_(f),REQ)is determined as a function of engine speed and accelerator pedal position from a lookup table stored in the ECU memory.

Given the estimation of intake composition from equation (32), the fuel supply can be limited to keep the resulting air/fuel ratio above a predetermined value to avoid the production of smoke. In step 108, this air/fuel ratio limit λ_(lim) is generated as a calibrated function of SOI(k), and one or more engine operating parameters such as engine speed, intake manifold pressure, engine coolant temperature, and fueling rate.

The exhaust burnt gas fraction (F_(e)), intake manifold burnt gas fraction (F₁), and predicted exhaust air/fuel ratio (λ_(e)), are given by: ##EQU4##

In this fuel control strategy, W_(c1) (k) is measured or estimated, W₂₁ (k) is calculated as shown in equations (12) (15), (26) or (30), and W_(1e) (k) is determined by equation (2). The requested fueling rate is also known. The variable T represents the number of time samples associated with the combined cycle delay and transport delay of the air/fuel mixture as it is inducted into the engine, combusted, exhausted into the exhaust manifold, and recirculated into the intake manifold through the EGR system. The time constant τ_(F21) represents the mixing of gases in the exhaust manifold and EGR path. Both of these time constants can be readily determined from engine mapping.

Alternatively, equation (34) can be substituted with a different filter given by the following equations:

    F.sub.1,pre (k)=F.sub.2 (k)W.sub.21 (k)/(W.sub.c1 (k)+W.sub.21 (k))(37)

    F.sub.1 (k+1)=F.sub.1 (k)+δt(1/τ.sub.F1 (-F.sub.1 (k)+F.sub.1,pre (k)))                                                     (38)

wherein τ_(F1) is a time constant representative of mixing of the gases in the intake manifold.

The performance of the estimated exhaust air/fuel ratio as described by equation (36) is shown in FIG. 3a. The simulated exhaust air/fuel ratio is shown by line 204 and the estimated exhaust air/fuel ratio is shown by line 206. It can be seen from the graph of FIG. 3a that the estimated exhaust air/fuel ratio tracked very closely with the simulated value with the exception that the simulation model became saturated at 90 and the estimator saturated at 120.

Step 112 compares the estimated exhaust air/fuel ratio and the air/fuel ratio limit necessary to avoid the production of smoke. If the predicted air/fuel ratio limit exceeds the limiting air/fuel ratio, the engine is operating in normal mode and the delivered fuel rate W_(f) (k) is set equal to the requested fueling rate (step 114).

If, however, the estimated exhaust air/fuel ratio is less than the air/fuel ratio limit, a change in the start of fuel injection (δSOI(k))timing is determined to be added to the standard injection timing (SOI(k)) to achieve a higher exhaust air/fuel ratio limit (λ_(lim),new (k)) (step 116). The change in injection timing is determined from a lookup table stored in the ECU memory and can be either a positive or negative value. The new exhaust air/fuel ratio limit (λ_(lim),new (k)) is then generated as a calibrated function of SOI(k)+δSOI(k), and one or more engine operating parameters.

The new air/fuel limit is then compared with the estimated exhaust air/fuel ratio in step 118. If the estimated exhaust air/fuel ratio is above the threshold represented by the new air/fuel ratio limit, the fuel injection timing is modified in step 120, and the fueling rate is set equal to the desired fueling rate in step 114. Otherwise the fuel is limited in step 122 according to the following equations:

    B=((1-F.sub.1 (k))W.sub.1e (k)/λ.sub.lim,new (k))-(F.sub.1 (k))W.sub.1e (k)/(Φ.sub.s +1))                        (39)

and

    W.sub.f (k)=min[W.sub.f,REQ (k),max(W.sub.f (k-1),B)]      (40)

Additionally, the ECU preferably commands the EGR system and VGT to increase the available fresh air intake to maintain the air/fuel ratio above the limit at which smoke occurs.

The performance of the novel fuel limited method is illustrated in FIGS. 4a through 4e. In FIG. 4a, line 400 shows the requested fueling rate in milligrams per stroke. FIGS. 4b and 4c show the resulting smoke in percent opacity (line 402) and exhaust air/fuel ratio (line 404) without the disclosed fuel limiting method. Line 406 of FIG. 4c represents the air/fuel ratio limit above which visible smoke occurs. FIGS. 4d and 4e illustrate the same variables at the same operating point where the novel fuel limiting scheme is used. As can be appreciated, when the fuel is limited (line 408 of FIG. 4a), the exhaust air/fuel ratio (line 410 of FIG. 4e) is maintained above the level at which visible smoke occurs (line 414 of FIG. 4e). Accordingly, the opacity of the smoke (line 412 of FIG. 4d) is maintained below the visible level.

While the invention has been described in connection with one or more embodiments, it will be understood that the invention is not limited to those embodiments. For example, although the engine system described includes a variable geometry turbocharger, the disclosed method would equally apply to engine systems with fixed geometry turbochargers or naturally aspirated engines as well. Accordingly, the invention covers all alternatives, modifications, and equivalents, as may be included within the spirit and scope of the appended claims. 

What is claimed is:
 1. A method of controlling the fueling rate of a compression ignition engine having an exhaust gas recirculation (EGR) system having an EGR valve connecting an intake manifold and exhaust manifold of the engine, and a turbocharger including a compressor and a turbine, the method comprising the steps of:generating a compressor airflow value (W_(c1) (k)) indicative of the mass airflow into the turbocharger compressor; generating a recirculation flow value (W₂₁ (k)) indicative of the flow of exhaust gas through the EGR system; generating an intake flow value (W_(1e) (k)) indicative of the flow from the intake manifold into the engine; generating an intake burnt gas fraction value (F₁ (k)); determining a start of fuel injection value (SOI(k)) as a function of engine speed and fueling rate; determining the requested fueling rate (W_(f),REQ) as a function of accelerator pedal position and engine speed; determining a predicted exhaust gas air/fuel ratio (λ_(e) (k)); generating an exhaust gas air/fuel ratio limit (λ_(lim) (k)); and delivering fuel to the engine (W_(f) (k)) as a function of λ_(e) (k), λ_(lim) (k), SOI(k), W_(f),REQ, F₁ (k) and W_(1e) (k).
 2. The method as set forth in claim 1 wherein the step of delivering fuel to the engine (W_(f) (k)) includes the step of comparing λ_(e) (k) to λ_(lim) (k) and, if λ_(e) (k)≧λ_(lim) (k) then delivering fuel to the engine (W_(f) (k)) at a rate equal to W_(f),REQ, else, generating a change in start of fuel injection value (δSOI(k)) resulting in a decreased exhaust gas air/fuel ratio limit (λ_(lim),new (k)) and, if λ_(e) (k)≧λ_(lim),new (k), then delivering fuel to the engine at a rate equal to W_(f),REQ, and modifying SOI(k) by an amount equal to δSOI(k), else, limiting the fuel supplied to the engine according to the following equation:

    W.sub.f (k)=min[W.sub.f,REQ (k),max(W.sub.f (k-1),B)]

wherein W_(f) (k-1) represents the previous fueling rate delivered to the engine and B is defined as:

    B=((1-F.sub.1 (k))W.sub.1e (k)/λ.sub.lim,new (k))-(F.sub.1 (k))W.sub.1e (k)/(Φ.sub.s +1))

wherein Φ_(s) represents the stoichiometric equivalence ratio.
 3. The method as set forth in claim 2 wherein the step of limiting the fuel supplied to the engine further includes the step of modifying the EGR flow and turbocharger flow to increase the airflow into the intake manifold of said engine.
 4. The method as set forth in claim 2 wherein the step of generating a compressor airflow value (W_(c1) (k)) includes the step of receiving an intake airflow signal from a mass airflow sensor.
 5. The method as set forth in claim 1 wherein the step of generating a recirculation flow value (W₂₁ (k)) includes the step of calculating the recirculation flow value as a function of the intake manifold pressure (p₁), exhaust manifold pressure (p₂), and EGR valve position.
 6. The method as set forth in claim 1 wherein the step of generating a recirculation flow value (W₂₁ (k)) includes the step of calculating the recirculation flow value in accordance with the following equation:

    W.sub.21 (k)=z(k)/T.sub.21 (k)

wherein z represents the scaled enthalpy flow estimate derived from an adiabatic assumption of the engine system and the first law of thermodynamics, and T₂₁ (k) represents the EGR system temperature.
 7. The method as set forth in claim 1 wherein the step of generating a recirculation flow value (W₂₁ (k)) includes the step of calculating the recirculation flow value in accordance with the following equation:

    W.sub.21 (k)=(1/T.sub.21 (k))((V.sub.1 /(γR))(dp.sub.1 (k)/dt)-f(k))

wherein T₂₁ (k) represents the EGR system temperature, V₁ represents the volume of the intake manifold, R represents the difference between the pressure specific heat constant and volume specific heat constant, γ represents the ratio of the pressure specific heat constant to volume specific heat constant, dp₁ (k)/dt represents an estimate of the time rate of change of the intake manifold pressure, and f(k) represents the filtered difference between the compressor mass air flow enthalpy and engine intake flow enthalpy.
 8. The method as set forth in claim 1 wherein the step of generating an intake burnt gas fraction value (F₁) as a function of W_(c1) (k) and W₂₁ (k) includes the step of determining the burnt gas fraction in the exhaust manifold (F₂); determining the intake manifold pressure (p₁); and calculating F₁ from a differential equation defined as:

    dF.sub.1 /dt=(RT.sub.1 /(p.sub.1 V.sub.1))[-(W.sub.c1 +W.sub.21)F.sub.1 +W.sub.21 F.sub.2 ]

wherein R equals the difference between the pressure specific heat constant and volume specific heat constant and V₁ represents the intake manifold volume.
 9. The method as set forth in claim 8 wherein the step of determining the predicted exhaust gas air/fuel ratio (λ_(e) (k)) includes calculating λ_(e) (k) in accordance with the following equation: ##EQU5## wherein Φ_(s) represents the stoichiometric equivalence ratio.
 10. A compression ignition engine system comprising:an engine block having a plurality of combustion chambers formed therein for combusting an air/fuel mixture; a plurality of fuel injectors corresponding to the plurality of combustion chambers, responsive to a fuel injection signal, for delivering fuel to said combustion chambers; an intake manifold for delivering intake air to the plurality of combustion chambers; an exhaust manifold for transmitting exhaust gas from the plurality of combustion chambers; a turbocharger having a compressor coupled to a turbine, the turbine being in communication with the exhaust gas in the exhaust manifold and the compressor being in communication with the intake manifold such that exhaust gas drives the turbine which causes the compressor to increase the flow of ambient air into the intake manifold; an exhaust gas recirculation (EGR) conduit having an EGR valve connecting the exhaust manifold and intake manifold of the engine for regulating the rate at which exhaust gas is recirculated into the intake manifold; a manifold pressure sensor in the intake manifold for providing an intake manifold pressure signal (p₁); a mass airflow sensor located upstream of the compressor for providing a compressor mass airflow signal (W_(c1) (k)); a temperature sensor in the intake manifold for providing an intake aircharge temperature signal (T_(c1)); and an engine control unit for generating said fuel injection signal comprising: a microprocessor programmed to:generate an intake airflow value (W_(1e) (k)) indicative of the airflow from the intake manifold into the combustion chambers; generate a recirculation flow value (W₂₁ (k)) indicative of the flow of exhaust gas through the exhaust gas recirculation system; determine a start of fuel injection value (SOI(k)) as a function of engine speed and fueling rate; determine the requested fueling rate (W_(f),REQ) as a function of accelerator pedal position and engine speed; generate an intake burnt gas fraction (F₁); determine the predicted exhaust gas air/fuel ratio (λ_(e) (k)); generate an exhaust gas air/fuel ratio limit (λ_(lim) (k); generate a fuel injection signal (W_(f) (k)) as a function of SOI(k), λ_(e) (k), λ_(lim) (k), W_(f),REQ (k), W_(1e) (k), and F₁ (k); transmit said fuel injection signal to said plurality of fuel injectors.
 11. The compression ignition engine system of claim 10 wherein the microprocessor generates a fuel injection signal (W_(f)) as a function of SOI(k), λ_(e) (k) , and λ_(lim) (k) by comparing λ_(e) (k) to λ_(lim) (k) and, if λ_(e) (k)≧λ_(lim) (k), then setting the fuel injection signal (W_(f) (k)) equal to the requested fueling rate (W_(f),REQ) else, generating a change in start of fuel injection value (δSOI(k)) resulting in a decreased exhaust gas air/fuel ratio limit (λ_(lim),new (k)) and, if λ_(e) (k)≧λ_(lim),new (k), then setting the fuel injection signal (W_(f) (k)) equal to the requested fueling rate W_(f),REQ, and modifying SOI(k) by an amount equal to δSOI(k), else, setting the fuel injection signal according to the following equation:

    W.sub.f (k)=min[W.sub.f,REQ (k),max(W.sub.f (k-1),B)]

wherein W_(f) (k-1) represents the previous fueling rate delivered to the engine and B is defined as:

    B=((1-F.sub.1 (k))W.sub.1e (k)/λ.sub.lim,new (k))-(F.sub.1 (k))W.sub.1e (k)/(Φ.sub.s +1))

wherein Φ_(s) represents the stoichiometric equivalence ratio.
 12. The compression ignition engine system of claim 10 wherein the microprocessor generates a recirculation flow value (W₂₁ (k)) by calculating the recirculation flow value in accordance with the following equation:

    W.sub.21 (k)=z(k)/T.sub.21 (k)

wherein z represents the scaled enthalpy flow estimate derived from an adiabatic assumption of the engine system and the first law of thermodynamics, and T₂₁ (k) represents the EGR system temperature.
 13. The compression ignition engine system of claim 10 wherein the microprocessor generates a recirculation flow value (W₂₁ (k)) by calculating the recirculation flow value in accordance with the following equation:

    W.sub.21 (k)=(1/T.sub.21 (k))((V.sub.1 /(γR))(dp.sub.1 (k)/dt)-f(k))

wherein T₂₁ (k) represents the temperature of the exhaust gas flowing from the exhaust manifold to the intake manifold, V₁ represents the volume of the intake manifold, R represents the difference between the pressure specific heat constant and volume specific heat constant, γ represents the ratio of the pressure specific heat constant to volume specific heat constant, dp₁ /dt(k) represents an estimate of the time rate of change of the intake manifold pressure, and f(k) represents the filtered difference between the compressor mass air flow enthalpy and engine intake flow enthalpy.
 14. The compression ignition engine system of claim 10 wherein the microprocessor generates an intake burnt gas fraction value (F₁) as a function of W_(c1) (k) and W₂₁ (k) by estimating the burnt gas fraction in the exhaust manifold (F₂), determining the intake manifold pressure (p₁); and calculating F₁ as a differential equation defined as:

    dF.sub.1 /dt=(RT.sub.1 /(p.sub.1 V.sub.1))[-(W.sub.c1 +W.sub.21)F.sub.1 +W.sub.21 F.sub.2 ]

wherein R equals the difference between the pressure specific heat constant and volume specific heat constant and V₁ represents the intake manifold volume.
 15. The compression ignition engine system of claim 14 wherein the microprocessor generates the predicted exhaust gas air/fuel ratio (λ_(e) (k)) by calculating λ_(e) (k) in accordance with the following equation: ##EQU6## wherein Φ_(s) represents the stoichiometric equivalence ratio. 